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ABSTRACT 

Energy release in the early Universe leads to spectral distortions of the cosmic microwave 
background (CMB) which in the future might allow probing different physical processes in 
the pre-recombination (z > 10 3 ) epoch. Depending on the energy injection history, the as- 
sociated distortion partially thermalizes due to the combined action of Compton scattering, 
double Compton scattering and Bremsstrahlung emission, a problem that in general is hard 
to solve. Various analytic approximations describing the resulting distortion exist, however, 
for small distortions and fixed background cosmology the Green's function of the problem 
can be computed numerically. Here we show that this approach gives very accurate results for 
general thermal histories, allowing fast and quasi-exact computation of the spectral distortion 
given the energy release rate. Our method can thus be used to forecast possible constraints on 
early-universe physics obtained from future measurements of the CMB spectrum. 
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1 INTRODUCTION 

The observations o f the cosmic microwave background (CMB) 
with COBE/FIRAS JMather et al.ll994l : lFixsen et alj 19961) showed 
that the CMB frequency spectrum is extremely close to a per- 
fect blackbody, allowing deviations of no more than A/ v //,, < 
few x 10~ 5 at wavelengths ranging from 100 /jm to 1cm. Dra- 
matic improvements in experimental design and technology over 
the past decades since the launch of COBE may soon allow much 
more se nsitive measurements of the CMB freque ncy spectrum 
(e.g., see iFixsen & Mathedl20C)2l ; iKogut et alj|201 ll) . For this rea- 
son, CMB spectra l distortion have recently experienced renewed 
interest (e.g.. see IChluba & Sunvaevl 1201 2|; [Chluba et alj|2012bl: 



Paier & Zaldarriaga]l2012l;lKhatri & Sunvaeyll2012bl;lTashiro et alj 
2012MAli-Ha mioud 2013l: ISunvaev & K hatri 2013, to name a few), 



potentially providing a way to probe new physics. It therefore is 
time to ask what exactly one might hope to learn about the ther- 
mal history of our Universe from future measurements of the CMB 
spectrum and how the information can be extracted. 

To address this problem, a fast and accurate method 
for computing the spectral distortion from different scenar- 
ios is required. Previou s numerical approaches (e.g., see 
IProcopio & Buriganal 120091 : IChluba & Sunvaevl 120121) are rather 
time -consuming. Simple analytic approximations exist (e.g 



see ISunvaev & Zeldovic hi Il970bl ; llllarionov & Sunvaev 



Burigana et alj 



1991 



k 



Hu&Silkl 1 19931 ; iBurigana et al 



1974; 



1995; 



Khatri & Sunvaev 2012b. a), however, for small distortions and 



fixed background cosmology it is possible to obtain a numerical 
approximation for the Green's function of the thermalization prob- 
lem. This allows us to compute the spectral distortion for general 
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thermal histories in an economic and quasi-exact way. In this short 
paper, we present our method which provides the basis for detailed 
forecasts, assessing the future experimental possibilities for spec- 
tral distortion measurem ents, and comparisons with independent 
JKhatri & SunvaewOlj) approaches. 



2 COSMOLOGICAL THERMALIZATION PROBLEM 

In the early Universe, photons undergo many interactions with 
the electrons and baryons. A detailed discussion of the cosmo- 
logical thermalization problem for differen t energ y-release mech- 
anisms can be found in lChluba & Sunvaevl J2012I) . Compton scat- 
tering allows photons to exchange energy with the electrons and 
baryons, causing diffusion in frequency, while double Compton and 
Bremsstrahlung emission are responsible for adjusting the photon 
number. The Boltzmann equation for the evolution of the average 
photon occupation number, n v (t), can be expressed as: 
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The second term on the left hand side describes the redshifting 
of photons due to the adiabatic expansion of the Universe and the 
right hand side terms correspond to the physical processes quoted 
above. We also explicitly added a photon source term, d/2 v /dr| s , 
however, direct photon production is normally negligible. If, for 
example, energy release is caused by some decaying or annihilat- 
ing relic particle, very few soft photons are directly produced, and 
hig h-energy photons and decay products just heat the medium (e.g. , 
see Shull & van Steenberg 1985; Chen & Kamionkowski 2004; 
iHiitsi e t al. 2009; Sl atver et alj 120091 : Icirelli et al] l2009h . Subse- 
quently, this leads to up-scattering of CMB photons and thus 



© 0000 RAS 



Chluba 



sources distortions, initiating the cosmological thermal ization pro- 
cess 

In addition to the photon occupation number, one has to follow 
the temperature of ordinary matter, including the cooling/heating 
due to the adiabatic expansion of the Universe, Compton scattering, 
double Compton emission and Bremsstrahlung. Extra energy re- 
lease can be parametrized using an effective heating rate, Q(t). It in- 
creases the matter temperature, subsequently causing up-scattering 
of CMB photons by Compton scattering and photo n production at 
very low frequencies ( ISunvaev & Zeldovichlll970al) . 

For us it is im portant that any pr i mordial spect ral distortion 
has to be tiny (cf.. iFixsen et al.lll996T . iFixsenl I20091) . so that the 
thermalization problem can be linearizeqj. We furthermore know 
the main cosmological parameters, such as the baryon density or 
Hubble expansion rate, extremely well from measurements of the 
CMB anisotropies JBennett et al.ll2003l ; |Planck Collaboration et al] 
12013ft . Thus it is possible to numerically compute the Green's func- 
tion, G t h(v, Zh, Zobs), of the thermalization problem using CosmoTh- 
erm. Here z h denotes the heating redshift and z obs the observing 
redshift. The resulting spectral distortion at z = is then given by 



AI v (z = 0) = I G,h(v, z , 0) — 



6zf 



(2) 



where p y (z) is the energy density of the CMB photon field, and 
we introduced the intensity, A/,, = (2hv 3 /c 2 )An v . We chose to use 
redshift as time variable, since most energy r elease mechanisms 
have rather simple form in this coordinate (see Ichluba & Sunvaevl 
I2012D . From the practical point of view this poses no significant 
limitation, and Eq. l|2} can be generalized to other variables, as long 
as the heating rate varies slowly. 



3 RESULTS FOR THE GREEN'S FUNCTION 

To obtain the Green's function we solve the thermalization problem 
assuming that the heating rate is given by a (^-function centered at 
heating redshift, Zh- Numerically we simply approximate the heat- 
ing rate by a narrow Gaussian in redshift 



d(Q/Pr) 



dz 



Ap y , h exp(-[z-z h ] 2 /2o" 2 ) 



Py.h 



(3) 



~ 6 . By evolving the distor- 
from z s » Zh until late 



witlyj <T C =a 0.01 and Ap y-h /p y .i, - 10 
tion, An v , and matter temperature, T, 
times (z — 200) we obtain an approximation for the Green's func- 
tion from Gth(v,Zh.0) as A«" um /[Ap rh /p r ,ii]. At z s the photon field 
is given by a blackbody at temperature T s < T (l + z s ), where 
T = 2.726 K jFixsedl2009h is the CMB temperature today. Ini- 
tially, also the temperature of ordinary matter is T m (z s ) = T s . At the 
end of the evolution, we ensure that the effective photon tempera- 
ture is close to To. To compute the Green's function we furthermore 
exclude the down-scattering of photons by adiabatically cooling 
electrons, leading to a distortio n characterized by negative chem- 
ical potential and y-par ameter dChlubal |2005t Ichluba & Sunvaevl 
120121 ; iKhatri et alj|20120 . This process can be included afterward, 
since for fixed cosmology the corresponding distortion is known 



1 We numerically solve the full non-linear system, however, limiting the 
total energy release to Ap y /p 7 <K 1. 

2 We confirmed the convergence of the Green's function by changing the 
width of the Gaussian, finding excellent agreement between the results. 



0.24, il m = 0.26, fi b = 0.044, fi A = 0.74, Q k =0,h= 0.71, and 
A'efr = 3.046, however, since the computation is rather straightfor- 
ward, it is very easy to adjust these values. 

3.1 Numerical results 

In Fig.[T]we present the numerical results for the Green's function 
of the cosmological thermalization problem. For very early energy 
release (zh » z M — 2 x 10 6 ), it is equivalent to a temperature shift, 
having the shape 

2hv 3 . 2hv 3 xe' 

4G th (v,z h » z„,0) * G(v) = -5- v<9,v/(v, T ) = — 

C c 1 (e l - \y 

where x = hv/kTo, n pl (v,To) = l/(e 1 ' - 1), and the factor of 4 
arises from Ap y jp y = 4AT/T. This indicates that the thermaliza- 
tion proces s is extremely efficient, with t he residual p-type spectral 
distortion ( ISunvaev & Zeldovichlll970ah being exponentially sup- 

S iressed by the distortion visibility functio r|j, ST(.Zh) = e - '- 1 '^''' 
Danese & de Zottill 19821 ; iHu & Silkll 19931) . To improve the accu- 
racy of our computation, at z > 5 x 10 4 we treated the temper- 
ature shift independently, imposing photon number conservation, 
fv 2 A« v dv = 0, to define the distortion parQ Since the tempera- 
ture shift caused by early energy release is not directly observable, 
this procedure allowed us to obtain precise results for the //-type 
distortion related to significant energy release at very early times. 

For low heating redshift (zh £ 1 4 ), the Green's function take s 
the shape of a Compton y-distortion l lZeldovich & Sunvaevlll969l) , 
At this redshift, the effective y-parameter remains very small, so 
that photons only weakly interact with the electrons and the distor- 
tion is described by 

4Gd,(v,Zh < 10 4 ,0) * Ksz(v) = —5- 7- .,, [xcoth(x/2) - 4] . 

c l (e J - \y 

At very low frequencies (v < 1 GHz) this approximation is 
not as accurate, since matter and radiation start decoupling and 
Bre msstrahlung emission (or absorption) can alter the spectrum 
(see Ichluba & Sunvaevll2012l for illustration). However, the total 
number of photons in that part of the radiation field is very small 
and this effect will be challenging to observe. The numerical results 
for the Green's function include these aspects. 

At Zh cs 3 x 10 5 the Green's function is mainly represented by a 
/j-type distortion, obtained using the condition f i^Any dv = (see 
IChluba et al]2012bHKhafri & Sunvaevll2012iJ . for more detail): 

2hv 3 e v 
a-'G lh (v,z h - 3 x 10 5 ,0) * M(v) = [x//3 - 1] . 

c 2 (e v - \y 



Here a « 1.40 1 JSunvaev & Zeldovicrj Il970d: 

llllarionov & Sunvaevl [l974l) and 3 = 3f(3)/f(2) ss 2.1923. 
Neglecting the temperature shift, the Green's function for 
the distortion part is well represented by M(v), even for 
Zh ^ 3 x 10 5 when additionally multiplying with the dis- 
tortion visibility function, jT(Zh). giving the approximation 



G th (y,z h > 3 x 10 5 ,0) * 1.401J r (zh)M(v) + 



i-J(--h) 
4 



G(v). 



3 See lKhatri & SunvaeM 120 12q) for recent analytical discussion. 

4 The thermalization process slowly increases the photon number density 
or equivalently the effective temperature (defined by N y oc T 3 ). To avoid 
having a large contribution from G(v) to the numerical solution we can 
restart the computation at different stages, shifting the temperature of the 
reference blackbody to the new effective temperature of the photon field. 
At z < 5 x 10 4 thermalization is already very inefficient and no significant 
average temperature shift can be produced, so that we leave the temperature 
of the reference blackbody unchanged. 
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Figure 1. Numerical results for the Green's function of the cosmological thermalization problem for various heating redshifts, Zh e [10 3 ,5 X 10 6 ]. Energy 
release at very high redshifts causes an increase in the effective temperature of the photon field, while at low redshifts photons only partially up-scatter, creating 
a y-distortion. Around Zh — 3 x 10 5 a pure ^-distortion is created. All intermediate stages are roughly (precision = 10% — 30%) represented by a superposition 
of these extreme cases, however, the small residuals provide a principle possibility to distinguish different thermal histories at redshifts 10 < z < 3 X 10 . 



Finally, Figure[T]also shows the Green's function at several in- 
termediate stages, 10 4 < Zh < 3xl0 5 . Si gnals produced mainly dur- 
ing this epoch were already discussed in lChluba & Sunvaevl (I2012L 
see Fig. 15 and 19), indicating that the total distortion is not sim- 
ply given by a pure supeiposition of f i and y-distortion. More re - 
cently, this was also demonstrated by iKhatri & Sunvaevl d2012ah . 
The small (^ 10% - 30%) residuals might allow distinguishing 
different energy release scenarios in the future. This is especially 
interesting if, for instance, decaying pa rticles with lifetimes f x - 
2.6 x 10 8 sec - 2.2 x 10" sec are present dChluba & Sunvaevll2012h . 
however, detailed forecasts are required to address this question. 
One can still obtain a fairly good approximation for the Green's 
function in this intermediate regime by assuming that the spec- 
trum is close to a superposition of a pure yu- and y-distortion, 
G th (v, Zh < 3 x 10 5 ,0) * 1.401 J^(z h )M(v) + ®2± Y sz (v). The 
coefficients, J y (Zh) and JT^(Zh), which approximate the transition 
between the /i- and y-distortion regime, can be determined using a 
least square fit to the Green's function. We find 



approximation for the Green's function at all redshifts is given by 



Jv(Zh) ~ 



1 + 



1+z 



6.0 x 10 4 



Jp(Zh) ~ 1 - exp 



1+z 



5.8 x 10 4 



(4) 



This approximation for G t h(v, Zh,0) works at =; 10% - 30% preci- 
sion for the standard cosmology, with the residuals containing addi- 
tional information about the time-dependence of the energy release 
at 10 4 < Zh S= 3 x 10 5 . In summary, this means that a pretty good 



G;(v,Zh,0) = 1.401J-„(z h ) JXzh)M(v) 



Jv(Zh) v , , 1-JXzh)^ , 
+ — -. — ^sz(v) + G(v). 



(5) 



The redshift-dependence of the functions, Sf y {lh)< J^(Zh) and J"(Zh) 
is illustrated in Fig.E] At z - 5.3 x 10 4 we find J y (z b ) - J^(zh), jus- 
tifying why z — 5 x 1 4 was previously used to mark the transition 
between u- and v-era JHu & Silk 1993). 

It is also important to mention that the energy integral 
/[G* h (v,zh,0)/ P) ,]dv = [j^(z h )J(zh) + j;v(zh) + l- l 7(zh)] < 1. 
The 'missing' energy can be found in the residual distortion, which 
is not represented well by the approximation, although it never ex- 
ceeds m 16% (maximum at z — 7.7 x 10 4 ). For simple estimates, 
Eq. lO should thus be very useful, while for high precision the ex- 
act form of the Green's function must be used. 



3.2 Spectral distortion for some examples 

We implemented the Green's function method as part of CosmoTh- 
eriv|j. With Eq. ((2) this allows very fast computation (a few sec- 
onds versus m 10 min - 1 h) of the resulting spectral distor- 
tion for given thermal history. In Fig. [2] we show the heating rate 
and resulting distortion for some illustrative cases. We confirmed 
the precision of the Green's function approach by explicitly com- 
puting the distortion for different scenarios solving the evolution 
equations with CosmoTherm, finding excellent agreement [^ 1%- 
level uncertainties in the background cosmology are more rele- 



CosmoTherm is available at www . Chluba . de/CosmoTherm 
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Figure 2. Effective heating rate and spectral distortion created by different scenarios. In all cases the total energy release was similar to Ap y /p y ^ few X 10~ 8 
(irrespective of observational constraints). We only included the distortion part for the calculation since the shift in the overall temperature of the photon field 
remains unobservable. For comparison we show a y-distortion with y = 10~ 8 . We also illustrate the redshift-dependence of J\.(zh). Ju(Zh) and S(Zh), which 
we multiplied by 10~ 9 , for convenience. The different scenarios are described in more detail in Sect. |3.2j 



vant]. The two annihilation scenarios assume different cross sec- 
tions, representing s-wave [(cry) a const], and p-wave ann ihi- 
lation [<ixv> =* T/m oc (1 + z)] (see iMcDonald et al J [20011. for 
related discussion). Strong constraints for s-wave annihilation of 
dark matter can be der ived from measurements of the CMB tem- 
perature anisotropies dGalli et all J2009I: iHiitsi et all l2009l 1201 ll ; 
IPlanck Collaboration et al.ll2013l) . however, CMB anisotropy lim- 
its on the p-wave scenario are much weaker because most of the 
energy is liberated prior to the recombination epoch (z — 10 3 ). We 



can see that for p-wave annihilation most of the energy is released 
during the /j-era, while for s-wave annihilation heating at low red- 
shifts also contributes significantly. Spectral distortions thus might 
allow shedding light on these cases, although degeneracies with 
other scenarios are expected. 

For illustration we also chose a decaying particle scenario with 
the main heating occurring at z - 5 x 10 4 . As explained above, 
in this case the distortion is not just given by a simple superpo- 
sition of (pure) /j and y-distortion, so that this scenarios might be 



© 0000 RAS, MNRAS 000, 000-000 



Green's function 



distinguishable from, for ex ample, the signal cr eated by the dis- 
sipation of acoustic modes JChluba et alj|2012b| ]a|), at least when 
the simplest shape for the primordial power spectrum is assumed. 
A more detailed forecast will be requireqj, and we plan to ap- 
ply our method, a ssuming experimental settings similar to PIXIE 
JKogutetalhoilh . 



3.3 Cosmology-dependence of the Green's function 

The cosmology-dependence of the Green's function is mainly 
driven by (i) the CMB monopole temperature; (ii) the number den- 
sity of baryons and electrons; (iii) the faction of helium atoms 
(though the Bremsstrahlung Gaunt-factors and the effect on the 
number of electrons per baryon); and (iv) the expansion rate of the 
Universe (i.e., N e e and at later stages fl m ). Since the cosmologi- 
cal parameters are known with percent-level precision, the Green's 
function is accurate to a similar level. For detailed parameter fore- 
casts the cosmology-dependence of the Green's function could thus 
lead to significant degeneracies, if differences at the level of a few 
percent matter. This can be easily taken into account, however, we 
leave an in depth discussion to future work. 



4 CONCLUSIONS 

We computed the Green's function for the cosmological fhermal- 
ization problem assuming standard cosmological parameters. With 
our method it is possible to predict the resulting spectral distortion 
signal for different energy release scenarios with high precision at 
practically no computational cost (see Fig [2] for some examples). 
This approach is now part of CosmoTherm and will be very use- 
ful for detailed forecasts of possible future constraints on the pre- 
recombination thermal history of our Universe derived from mea- 
surements of the CMB spectrum. 

In particular, it will be very important to ask how degen- 
erate different scenarios are. For example, energy release related 
to dark matter annihilation with s-wave cross section is expected 
to give rise to a distortion that is practically degenerate with 
those caused by dissipation of adiabatic modes in the early Uni- 
verse, if the spectral index n s of the primordial power spec- 
trum at wavenumber k > lMpc~' is scale-invariant (e.g. see 
IChluba & Sunvaevl 20121). One the other hand , recent constraints 
from Planck ( IPlanck Collaboration et al .11201 31) show that at large 
scales n s =* 0.96. Extrapolating this to small scales thus suggests 
that the shape of the distortion ought be different in the two cases. 
Our computations also show that thermal histories with decay- 
ing relics might be constrainable if the lifetime of the particle is 
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t x — 2.6 x 10 s sec - 2.2 x 10" sec. We look forward to exploring 
these possibilities in the future. 
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